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ABSTRACT 

We present the first strong lensing study of tlie mass distribution in tlie cluster 
MS 2053-04 based on HST archive data. This massive, X-ray luminous cluster has 
a redshift z=0.583, and it is composed of two structures that are gravitationally 
bound to each other. The cluster has one multiply imaged system constituted 
by a double gravitational arc. 

We have performed a parametric strong lensing mass reconstruction using 
NFW density profiles to model the cluster potential. We also included pertur- 
bations from 23 galaxies, modeled like elliptical singular isothermal sphere, that 
are approximately within V x 1' around the cluster center. These galaxies were 
constrained in both the geometric and dynamical parameters with observational 
data. Our analysis predicts a third image which is slightly demagnified. We 
found a candidate for this counter-image near the expected position and with 
the same F702W-F814W colors as the gravitational arcs in the cluster. The re- 
sults from the strong lensing model shows the complex structure in this cluster, 
the asymmetry and the elongation in the mass distribution, and are consistent 
with previous spectrophotometric results that indicate that the cluster has a bi- 
modal mass distribution. Finally, the derived mass profile was used to estimate 
the mass within the arcs and for comparison with X-ray estimates. 



-2- 



Subject headings: Gravitational lensing: strong lensing — Galaxies: clusters — 
Galaxies: clusters: individual(MS 2053-04) 



INTRODUCTION 



The possibility that clust ers of galaxies might produce gravitational lensing had been 



1984I ). Twenty years ago (iLynds fc Petrosian 



consi dered by several aut hors (IKaroii fc Nottalelll976l : iDyer fc Roedei 



19861 ) and (ISoucail et al. 



1976 



Narayan et al. 



19871 ) independently 



announced the detection of an arc-like structure around the center of Abell 370, but it was 



Paczyhskil (119871 ) who explained the phenomenon in terms of gravitational lensing. Currently, 
the strong and weak regimes of the gravitational lensing effect are often used for constraining 
the density profiles of galaxy clusters. 

The massive clusters of galaxies have a surface density la rger than the critical value for 
strong lensing, approximately 1.0 g cm~^ (ITurner et al.lll984l ). They can produce multiple, 
strongly elongated images (giants arcs) of background galaxies lying on top of the caustic 
lines generated by the cluster potential. Observati onally, the first dedicated search for grav- 
itational arcs in X-ray clusters was performed by iLe Fevre et al.l (Il994j ) using the Einstein 
Observato ry Extended M e dium Sensitivity Survey (EMSS). A subsequent search was per- 
formed by iLuppino et al.l (119991 ) . which allowed to confirm the relationship between X-ray 
luminosities and cluster masses; these authors found strong lensing in eight of 38 clusters. A 
subsequent survey based on optically selecte d clusters in the Las Camp anas Distant Cluster 



Survey confirmed these giant-arc fractions (IZaritsky fc Gonzalea l2003l ) and iGladders et al. 



(120031 ) reports even larger fractions in the Red-Sequence Cluster Survey, showing that giant 
arcs are quite common in clusters of galaxies. 

Since the shapes of these arcs are related to the gravitational potential, the strong 
lensing offers a method to investigate how the barionic and dark matters are distributed 
in the inner regions of the cluster. Also it is possible to put significant constraints on 
the cosmological parameters using a set of stron g lensing clusters showing multiple images 
with spectroscopic redshifts (ISoucail et al.ll2004l ). Mass distributions have been modeled for 
many clusters with giants arcs. In a few cases, accurate models for the mass distribu tion 
have led to the identification of additional counter-image candidates (IKneib et al.lll996l ). In 
other cases, the detectio n of the predicted i r nages is uncertai n, since they are lost in the 
cD light distribution (see lCavazzi et al.ll2003l : ISand et al.l 12004 ). Recently, the gravi t ationa l 
lens cluster A 1 6 89 ha s been modeled in unprecedented detail by iBroadhurst et al.l (120051 ): 
Limousin et al.l (120061 ). These authors combined HST Deep Advanced Camera images and 
extensive ground based spectroscopy, identifying 32 multiply lensed systems. 
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The aim of this work is to describe the density profile and the mass distribution of 
MS 2053, using strong lensing. This cluster has a redshi ft of z=0.583, is not very optically 
rich in comparison with other z > 0.5 EMSS clusters (iLuppino et al.lll999l ). and among 
these high-redshift clusters it has the lowest X-ray luminosity. BeppoSAX observation s 



yield a X-ray luminosity Lx(2-10 keV) = f3.9±1. 0)xl0^^ erg s-^/lDella Ceca et al.ll2000f) 



and C handra (0.5-2 keV) = 3.5 xlO^^ erg s"^ Jvikhlinin et al.ll2002h . Ihuppino fc Gioia 
(I1992I ) discovered a gravitationally lensed blue arc in deep images of MS 2053, located ~ 16" 
in the North direction from the Brightest Cluster G alaxy (BCG). The arc is ~ 11" long and 
breaks into two clumps, labeled A and B. Recently, iTran et al.l ( 120051 ) were able to measure 
the redshift of image A using Keck spectroscopy data and determining the source redshift 
to be z = 3.146. 



Hoekstra et al.l (120021 ) detected the weak lensing (WL) signal induced by this cluster 
of galaxies from a two-colour mosaic of six Hubble Space Telescope (HST) Wide Field and 
Planetary Camera 2 (WFPC2) images. They fitted a singular isothermal sphere model to 
the observed azimuthally aver aged tangential distor sion and estimated the cluster velocity 



dispersion a = 886 km s ^. iHoekstra et al.l (l2002l ) also fitted the predicted profile from 



the NEW halo to the observed tangential distorsion of MS 2053, and found a concentration 
parameter (see section 3.1) in good agreeme nt with the pred icted value for a cold dark 
matter model without cosmo logical constant. iTran et al.l (120051 ) combined images from the 
wide-field HST WFPC2 and extensive spectral data from the Keck Low Resolution Imager 
Spectrograph (LRIS) to make a detailed study of the galaxy populations in MS 2053. They 
found that the cluster is composed of a main cluster (MS 2053-A) and an infalling structure 
(MS 2053-B) that are gravitationally bound to each other; their study also shows that the 
extended spatial distribution and lack of spectral segregation in MS 2053-B is consistent 
with a scenario in which it is in the initial stages of being accreted by the main cluster. Of 
the 149 spectroscopically confirmed cluster members, 113 belong to MS 2053-A and 36 to 
MS 2053-B. From these previous researches we present in Table [T] a compilation of physical 
quantities relevant to our work . 



Despite all the research cited above, only iHoekstra et al.l (|2002| ) have made an attempt 
to describe the strong lensing effect for this cluster and it was done when the arc redshift 
had not been measured yet. Assuming a redshift of z=2 for the arc and adopting a singular 
isothermal sphere, these authors found a velocity dispersion of about 1030 km s~^, consistent 
with their WL estimate. They argue that the mass distribution is elongated in the direction 
of the arc, which can be very well appreciated in the light distribution, and this elongation 
is the cause of the high value of their strong lensing mas s estimate. Such e longation in the 
distribution of the cluster members was also reported by lTran et al.l (120051 ). To assess with 
precision this mass distribution, we present a detailed strong lensing model of the bimodal 
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cluster MS 2053. 

In §[2|, we describe the observational data and present our sample of field galaxies. In 
§ini we define the model of the lensing cluster, and depict the profiles for both the cluster 
and the Galaxy-scale mass components. In §IH we present the main results. In §[5l we 
summarize and discuss our results. Finally, in §[6l we present the conclusions. Throughout 
this paper, we adopt a spatially-flat cosmological model dominated by cold dark matter and 
a cosmological constant. We use flm = 0.3, = 0.7 and Hq = 70 km s~^ Mpc~^. With 
these cosmological parameters and at the redshift of MS 2053, 1" corresponds to 4.62 
kpc. 



2. THE OBSERVATIONAL DATA 



The HST data have been obtained from the Multimission Archive at Space Telescope 
(MAST). They consist of two WFPC2 imagetQ obtained in both the F814W and the F702W 
filters; the total integration time for each filter was 2600 and 2400 s, respectively. The 
image reduction was performed using the IRAF/STSDASS package. First, a warm-pixel 
rejection was applied to the images using the IRAF task warmpix. The cleaned images were 
then combined with the task crrej to remove cosmic-rays hits. Finally, the background was 
subtracted and the WFPC2 chips were combined using the task wmosaic. In Figure [T] we 
show a 60" X 60" field obtained from the final four-chip mosaic frame for the F814W filter; 
this is the region used in our analysis of strong lensing. 

The magnitudes and geome tric parameters of the galaxies were measured using SEx- 
tractor (iBertin fc Arnoutd Il996l ). For the arc photometry we used a different p rocedure 
because SExtractor often overestimate the sky background (see lSmith et al.ll2005l ). and the 
arcs are highly distorted. These two circumstances reduce the precision of automatic pho- 
tometric measurements. Therefore, we employed polygonal apertures to obtain more accu- 
rate measurements. The vertices of the polygons for each arc were determined using the 
IRAF task polymark, and the magnitudes inside this apertures were calculated using the 
IRAF task polyphot. Hence, we found F814W magnitudes of mpgi^w = 21.64 ± 0.02 and 
^F8MW = 22.01±0.04 for arcs A and B respectively, and therefore the total magnitude for the 
system is 21.06 ± 0.11. Similarly, the F702W magnitudes are mp7Q2W = 22.05 ± 0.02 for arc 
A and TaFjQ2w = 22.38 ± 0.04 for arc B, contributing to a total magnitude of 21.45 ± 0.11. 
Our measurements of the tota l brigh tness in these two bands are slightly lower than the 
values reported by ISand et al.l (120051 ) (jaFsuw = 20.90 and inFio2W = 21.27). The small 



^Proposal ID 5991 



- 5 - 



discrepancies probably ari se from the differ ent metliodologies used to calculate the total 
magnitudes. For example, ISand et al.l (120051 ) used a single aperture to measure the magni- 
tude for the system of arcs as a whole (arc AB), while we used two apertures to measure 
each arc independently. 

We also calculated the F702W-F814W color for the AB arc system. The color for arcs 
A (0.41 ±0.03) and B (0.37±0.06), are similar within the errors, ar id the color for the whole 
AB system, 0.39 ± 0.15, is similar to the value of 0.37 reported by ISand et al.l (120051 ). 



The general properties of the cluster members are pre sented in Ta ble [2l We used 
the relative positions for the individual galaxies provided by iTran I (120021 ) to calculate the 
coordinates of each member of the cluster. Column (1) lists the identification for each 
MS 2053 galaxy. Columns (2) and (3) list the right ascension and declination. Column (4), 
the F814W magnitudes. Columns (5) and (6) show the geometric parameters derived from 
SExtractor: ellipticity e and the position angle ^o, which gives the orientation of the semi- 
major-axis from the horizontal line in the image, measured counter-clockwise. Column (7) 
lists the redshifts. Finally, columns (8) and (9) list the value s and the referen ces for the 
central velocity dispersions. The upper limits were taken from iTran et al.l (120051 ) (see their 
Figure 11), who consider only galaxies brighter than M^e < —20 + 5logh. 



3. THE MODEL OF THE LENS 



3.1. Dark matter component 

For studying the mass distribution in MS 2053, we use standa rd lens modeling techniques 
implemented in the LENSTOOlfl ray-tracing code (jKneiblll993l ). 



We model th e cluster with a NFW profile that has been predicted in cosmological N- 
body simulations (INavarro et al.l 119971 ): 



p(r) 



Pcrit{z)Sc{z) 

{r/rs){l + r/r-s 



where 5c is a characteristic density contrast, Pcrit{z) is the critical density of the universe as 
a function of the redshift z, and is a scale radius that corresponds to the region where 
the logarithmic slope of the density equals the isothermal value. However, higher-resolution 
simulations have suggested that the inner cusp of the NFW profiles is too shallow (e.g. 



^ This software is publicly available at: |littp://www.oamp.fr/cosinology/ 



lenstool/ 
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Moore et al.l Il998l : iGhigna et a. 



It is worth noting that iGavazzi 



20001) . giving origin to the generahzed NFW-type profile. 
2005h studied MS 2137 - 2353, obtaining similar results by 



using both profiles for the strong lensing regime. Therefore, for simplicity we adopt a profile 
similar to the original NFW. However, we consider an elliptical profile since in modeling the 
gravitational lens it is important to take into account the asphericity in the cluster potential. 

To generalize the spherical ra o del to an ellipsoidal one, we used the "pseudo-elliptical" 
NFW proposed by lGolse fc KneibI (120021 ). This potential is characterized by six parameters: 
the position X,Y; the ellipticity e; the position angle 9; and the velocity dispersion ag. 
We define ag as: 



(J. 



4 1 + In 



Gr'^gpcritSc 



(2) 



This is approximately half of the value proposed by iGolse fc KneibI (120021 ) and it is defined 
in such a way that it represents a realistic velocity dispersion at radius and not only a 
scaling parameter (see the appendix). 

We can relate this ag to the virial mass Ma and the virial radius ta of the cluster. If 
we define ta as the radius of a spherical volume within which the mean density is A times 



the critical density at the given redshift z (Ma 
expressed by: 



^77ir\pcrit^) , then the virial radius can be 



rA = [2MAG/AH{zf] 



211/3 



(3) 



The value of the characteristic overdensity A is taken from the solution to the collapse of a 
spherical top-hat m odel. Its value depends on fl and can be ap proximated by A = nSQ'^^ 
when Qjn + = 1 ( ILacey fc Cold 1 19931 : lEke et al.l 1 1999 . Il998l ). Integrating equation [1] and 
using the above definition, it is straightforward to show that the ratio between virial and 
scale radii, which is commonly called the concentration, c = r^/rg, is related to 6c by 



3 + c) - T^] 



(4) 



Solving equation m we can determine the concentration value given the z, ag and r^. 
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3.2. Galaxy-scale mass components 



For simplicity, we modeled the 23 galaxies using elliptical si ngular isothermal spheres 
(ESlSs) for the individual galaxies (IBlandford &: Kochanekl 119871 ) . This galaxy-scale mass 
components were included as perturbation into the cluster potential. The potential for ESISs 
is given by 



i^{r,e) = Arc 



D 



LS 



c2 



OS 



Vl + ecos(2(^-0o)) 



(5) 



where ctq corresponds exactly to the true 3D velocity dispersion and 6*0 is the position angle. 
Note that equation [5] is defined completely if we can obtain the values for ctq, the ellipticity 
e, and from observational data. The luminosity distribution of a given galaxy does not 
trace the dark mass distribution in its halo, however there is ev id ence that the projected 
mass and light distributions tend to be aligned (IKeeton et al.lll998l ). iKoopmans et al.l (120061 ) 
report this alignment between light and mass for a sample of early-type lens galaxies, for 
which the isophotal and isodensity countour trace each other. These authors also report 
an alignment between the stellar component and the singular isothermal ellipsoids used in 
their lens models. Following these authors, we infer the parameters of the galaxy-scale 
mass components from the position, the ellipticity and the position angle of the luminous 
component. 



Wuyts et al.l (120041 ) have r neasured the veloc ity dispersions of several elliptical and lentic 



ular galaxies in the cluster, and lTran et al.l (120051 ) have estimated the internal velocity disper- 
sion for the whole cluster sample with M^e < —18-1-5 log /i. Nevertheless, the individual data 
for each galaxy has not been published for all the objects. To minimize the number of model 
free parameters, the galaxies for which the velocity dispersion is not avail able (see Table [2|) 
are s caled as a function of luminosity following the Faber- Jackson relation (IFaber fc Jackson 
1976h : 



^0 = (Tq^ — 



1 1/4 



(6) 



which transformed to magnitudes is expressed by loguo = —0.1m + K, where K is a scaling 
factor. In this way the cluster members are incorporated into the model as a galaxy scale 
perturbation with all their individual parameters fixed. 
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3.3. The model optimization 

In principle, the optimizations in the source plane and in the image plane are mathe- 
matically equivalent. However, we chose using the source plane because the solutions are 
numerically simpler and faster to compute. The fit of the model is optimized by mapping 
the positions of the multiple images back to the source plane and requiring them to have 
a minimal scatter. To achieve this goal we use a figure-of-merit-function, x^, to quantify 
the goodness of the fit for each trial of the lens model i.e., a function that measures the 
agreement between our data and the fitting model for a particular choice of the parameters: 



2 _ 2 _|_ 2 /y\ 



The first term is constructed as follows: given a model, we compute the position in the 
source plane uf = {xf,yf) for each observed image {I < i < N), using the lens equation. 



— '^^ — Vip{ul), where is the lens potential. Therefore, the barycenter = {x^,y^) is 



constructed from the N sources, and the x for the position is defined by: 



2 _ 1 ''^ {x^-xf)'^{y^-yT)- 

^pos fJ 2-^ ^2 



were a^os is the error in determining the position of the source. 

The second term, xihape^ measures how well the shape of the arcs is reproduced by 
the model. Given a model, we compute the complex deformation, fj, for each observed 
image (1 < i < iV ), using the lens equation for complex deformation (see the discussion in 



Kneib et al.l 119961 . and their detailed derivation). The complex deformation of the image 
represents how much the intrinsic source shape changes due to the induced deformation 
caused by the potential. This complex deformation is function of the ellipticity, the position 
angle and the lens potential, f/ = ^/(ef , 9^ , ip). It has an equivalent form in the source plane, 
ff = ff {ej , 9j , (f) . If we define r^j and r^j as the components of the complex deformation 
in the source plane, and aghape as the error in determining the shape of the source, we get: 

2 ^ _ l^x',i T^x,i-l) + yy,i ^y,i-l) 

Xshape Ty / ^ 2 \ ) 

j^—l shape 

The errors, apos and a shape, have to be measured in the source plane, but these are not directly 
observable quantities. To avoid this problem,we assumed that the errors in the source plane 
arise from the measurements performed in the image plane. Thus, we arbitrarily set apos and 
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o'shape to 0.2" and 0.8" respectively, because these are the accuracies achieved in the position 
and shape measurements from the HST images. This is quahtatively correct in the sense 
that for poor models, the merit function (equation [7]) gives large values and thus indicates a 
poor fit. 

In the 6.2 version of LENSTOOL the Bayesian Markov Chain Monte Carlo (MCMC) 
package BAYESI^ has been implemented. The Bayesian MCMC method will be described 
in a dedicated forthcoming publication (Julio et al. 2007, in prep.) This package allows to 
obtain the errors and to avoid local minima. According to a user defined model parameter 
space, the MCMC sampler draws random models and compute their x^- Progressively, it 
converges to the most likely parameter space and outputs the best solution. We used these 
random models to perform the statistics to compute error bars on the estimation of the 
parameters, and to constrain the positions of the bimodal cluster (see next section). 



4. RESULTS 
4.1. A simple model for arc A 



The only arc for which the redshift has been measured is arc A fiTran et al.ll2005l ). The 
first step was to demonstrate that arcs A and B belong to the same source. But a model that 
includes only arc A and not arc B is unconstrained, and the data available is not adequate 
to break the degeneracy. Given the general characteristics of the cluster, and a gravitational 
lensed source with the properties inferred from arc A (i.e., redshift and arc position), we 
investigate if this source would produce multiple images or not. To achieve this task, we 
constructed a simple spherical model with a NFW profile and its central 23 galaxies modeled 
like elliptical isothermal spheres. 

Using equation H] we calculated the virial radius of the cluster. For a virial mass of 
3.7 X lO^^M© and a cluster redshift of z=0.583 (see Table [1]), ta ~ 2.9 Mpc. Therefore, 
~ 1.0 Mpc with a concentration parameter c=2.9 (according to the prescription described 



in section 3.3 of lEke et al.ll200ll ). The characteristic velocity, cXs = 1600 km s~^, is obtained 
from equation [2J With these values (r^ and as), the NFW profile is fully characterized. At 
this stage, the ellipticity and the position angle are set to zero. 

The velocity dispersion in the galaxy-scale mass components was fixed using the scaling 
relation given by equation [61 We found that the early-type cluster galaxies in the sample of 



3j. Skilling, 



http : / /www. inference, phy. cam. ac .uk /bayesys / 
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Wuyts et al.l (12004 ) satisfy the Faber- Jackson relation with log an = — (0.1± .05)m + (4.4 ± 
1.1). For the remaining subset of galaxies not included in the IWuyts et al.l (120041 ) sample, 
only galaxy n has a measured (Tq value. We supposed that these galaxies also satisfy the 
relation but with a different scaling factor K. In order to construct the new correlation, 
we considered the luminosity L* corresponding to the faintest galaxy in Table [2] (galaxy w , 
niF8i4iy = 22.57). We set = 64.7 km s~^ such that the galaxy n satisfies this correlation. 
With these v alues , the velocity dispersion in the subset of galaxies not included in the 
Wuyts et al.l (120041 ) sample is expressed by logdo = —0.1m + 4.19. Figure [2] sh ow both, the 



previo us logao and the Faber- Jackson relations for the early- type galaxies of IWuyts et al. 



With all these parameters we calculate the counter-images of arc A. The model obtained 
with these fixed parameters predicts two more images, one very demagnified and the other 
only one magnitude bellow the arc A; the position does not agree with B, but this shows that 
for this cluster characteristics and with this redshift of the source, A would have multiple 
counter-images. Therefore, is entirely plausible that images A and B arise from the same 
source. 



4.2. Bimodality evidence 

Two kinds of models are possible depending on whether we as sume that MS 2053 is 
composed of one or two clusters, as reported by iTran et al.l (120051 ). We performed the 
fits setting up the position (two parameters), the velocity, scale radius, ellipticity, and the 
position angle of the cluster as free parameters. These 6 parameters define the lens potential 
and they comprise all the allowed free parameters for this single cluster case. With just 2 
images and 4 constrains (if there is a total of Yl^=i = ^ images, then there are ^"^^ 4(?T,j — 
1) = Nc constrains on the models assuming that the position, ellipticity, and position angle 
of the images are fitted) the model is underconstrained. And it gets worse in the second 
model where we have 12 parameters defining the lens potential. Therefore, we expect that 
there may be a family of models that can fit the data, even with unrealistic parameters. 
Trying to avoid this problem, we have limited one parameter in the NFW, the scale radius 



r.«. 



The Tg cut-off was made in the following way: Given the cluster redshift and using 
equation [3|, we can write the scale radius as r<j oc being M the mass found from 

lensing. Since lensing is sensitive to the integrated mass along the line of sight, we expect 
discrepancies bet ween the t rue 30 mass, Mtme, and the lensing mass estimates of our models. 
Miens- Following iGavazzil (120051 ) . we consider an axisymmetric NFW density profile with 
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either the major or the minor axis ahgned toward the hne of sight. Therefore, the net effect 
of the projection is given by Mi^ns = QMtrue] where q is the ratio between the semiaxis 
ahgned and the semiaxis perpendicular to the hne of sight. Thus, we can write: 



(10) 



Using the concentrations shown in Table [T] and the axis values from iGavazzil (120051 ) we can 
constrain the scale radius for each model. For the one component model, 0.60 Mpc < < 
1.74 Mpc. Similarly for the two component model 0.30 Mpc < < 0.87 Mpc and 0.10 Mpc 
< < 0.28 Mpc for component one and two respectively. 



We constructed three models for the gravitational lens system. The first model, which 
is unconstrained, consists of a single cluster and data from arcs A and B to make the 
fit. Our best fit for the single component model predicts a third image which is slightly 
demagnified {mpsuw ~ 23) in a position where no object is detected in the HST exposures. 
We explored the field near its predicted position and identified an object around 6" to the 
north (R.A. 20° 56' 22'.'5 and DEC. - 04^ 37°^ 37*^) which seems to be a deformed image of 
a galaxy. Fig. [3] shows the isocontour map superimposed onto the optical image of this 
blurry object. We measured the magnitude of this counter-image candidate (jtifsuw = 
23.37 ± 0.05), and we found that it is similar to the magnitude predicted by our model. 
Moreover, its F702W-F814W color turns out to be analogous to arc A, 0.46 ± 0.07, (see 
discussion below). 

In the second model, we repeated the previous fit but including the third image candi- 
date. With this third image the model is well constrained because we have 8 constrictions 
and only 6 parameters. The fit was very accurate for the position but the total is affected 
by the shape. It is interesting to notice that in this well constrained model, there are two 
different solutions in the XY-space. In Figure H] we show the XY parameter space for this lens 
model constructed from all the trials coming out from the MCMC optimization procedure. 
We think that this result is due to the fact that the lens is a double cluster; we used the 
regions defined by these clumps to construct a third model, the double component model, 
that accounts for a bimodal distribution of mass. 

The results of these fits are summarized in Table [31 Column 1 identifies the model: 
SINGLE for models with one component and DOUBLE for the model with two components; 
AB and ABC indicate that the model consists of arcs A and B, or these two arcs plus the third 
image C. Columns 2 and 3 show the position in arcseconds relative to the BCG. Columns 4-8 
list the parameters associated to NFW profile, ehipticity, position angle, concentration, scale 



radius and a^, respectively. Column 9 and Column 10 show the ^ot the image positions 
only and for the whole fit, respectively. 



4.3. Mass profile 



For the purpose of quantifying the mass distribution, we traced a projected mass map 
from the best-fit model. Double ABC (Figure [5]). Figure [6] shows the corresponding iso- 
contour map superimposed onto the F814W image of MS 2053. By integrating this two 
dimensional map, we can get the total mass of the cluster. In other words, we can determine 
the mass, making an azimuthal average of its mass in circles around the BCG. 

In the figure [7| we compare the projected mass as a function of radius for the three 
components of the cluster: The two clusters {Component 1 and Component 2) and the 
individual galaxies. To calculate the contribution of the individual galaxies, we extrapolate 
the masses of the individual cluster components to their virial radii and add them. We 
can appreciate that within the central 0.8 Mpc, the dominant component corresponds to the 
mass of the galaxies. However, at the arc radius (r ~ 0.1 Mpc) the Component 1 has the 
main contribution to the lens mass. Figure [7] also show the total mass of the cluster as the 
sum of these three components. 

Figure [8] shows the sum of the projected mass of Component 1 and Component 2 (with- 
out including the galaxies) as a function of radius. We have extrap olated the NFW p rofiles 
up to 1 Mpc in order to co mpare with the virial mass estimate by iTran et al.l (120051 ). and 
with a j3 model derived by IVikhlinin et al.l (120021 ). The vertical line shows the distance of 
arc A with respect to the BCG, and it indicates the limit below which the mass estimate 
from strong lensing is reliable. The mass determined using X-ray data depends on assump- 
tions involving spherical symmetry and hydrostatic equil ibrium. Howeve r, MS 2053 shows 
evidence of being two clusters enduring a merger process (ITran et al.ll2005l ). and therefore it 
is not expected to find a good agreement between the (3 model and the strong lensing mass 
profile (see §|5]). 



5. DISCUSSION 

We found that the colors for the ABC arc system turned out to be virtually the same, 
with arc B slightly bluer but still coinciding with the color of arc A and image C within the 
errors. This result supports the assumption that images A, B and C belong to the same 
source. Despite this fact a serious concern, is the lack of specularity; arcs A and B should 
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be mirror images of each other, for which there is no evidence in the HST image. However, 
this would not be the first case of lack of specularity. For example, the cluster A 383 shows a 
giant tangential arc with substructures and with no evidence of being formed by two mirror 
images; these sub structures have d i fferent colors (ISmith et al.l I2OOII ) and he at the same 
redshift, z = 1.01 JSand et al.ll2"ooi ). ISmith et all J2005h have shown that these arc features 
can be understood in terms of just one galaxy positioned over the tangential and radial 
caustic lines simultaneously, with three different portions of the same galaxy producing the 
arc substructure. We think that something similar may be occurring in MS 2053, that is the 
proximity of the source to the caustic line (see Figure [9]). The color distribution in galaxies 
often changes at different optical radii (e.g., bulge, disk and halo populations). In the case 
of a spiral galaxy for which the disk (but not the bulge) is located over a caustic line, part 
of the disk is mapped onto a stretched arc such as arc B. This effect may explain the slight 
difference, if any, for the arc B color, the difference in shape between arcs A and B, and the 
lack of specularity between the arcs. 

For the rest of this section we will assume that the three images A, B, and C correspond 
to the same source. This assumption, along with the evidence presented in the following lines, 
in the sense that the quality of the fits is better for the model with two cluster components, 
justifies that our discussion shall be focused on our model 3. The for the whole fit is < 2.0 
for the DOUBLE ABC model and > 3.0 for the other two, favoring the bimodal model (cf. Xtot 
values in Table [3]). In addition the as values, 746 ± 114, and 392 ± 133 for Comp onent 1 and 



Comp onent 2 respectively, agree within the errors with the values reported by iTran et al. 
( I2OO5I ): but the other models yield larger velocities {cf. Tables [T] and [3]) . 



Hoekstra et al.l (120021 ) suggested that the mass distribution in MS 2053 is elongated in 



th e direction o f the g iant arc; this projected spatial distribution of the cluster was confirmed 
by iTran et al.l (120051 ) and it was attributed to a real elongation in the spreading of cluster 
members. The critical lines for our model, shown in Figure [HI reveal that the main component 
of the double cluster has a position angle of ~ 80°, and supports the existence of a mass 
elongation. In addition, the mass map as well as the isocontour map of the projected mass 
distribution (see Figure [5] and Figure [6]) refiect MS2053's elongated distribution. The offset 
in the position of the cluster components with respect to the BCG (see Table [3]) also confirms 
the c luster asymmetry. Similar offsets have been found by other authors for different clusters 



(.e.g. 

in unrelaxed clusters ( 



Smith et al. 



200 



Oguri et al.l 12004 : ICovone et al.l 120051) . and the effects can reach values up to 120 kpc 



Several researches hav e reported large concentration values by fitting NFW profile s 
to clusters of galaxies (e.g. iGavazzi et al.l l2003l : iKneib et al.ll2003l : iBroadhurst et al.ll2005l ). 
However, these results were obtained assuming spherically symmetric mass distributions, 
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althoug h the CDM model predicts triaxial halos as co nsequence of the c ollisionless dark 
matter (jjing &: Sutd l2002l ) . Using N-body simulations, IClowe et al. J2004h argue that the 
halo triaxiality af fects the conce ntration parameter measurements obtained through grav- 
itational lensing. iGavazzil (120051 ) concludes that a prolate halo aligned toward the line of 
sight is a natural explanation for the high concentration found in MS 2137-2353 ; but he also 



l eaves open the possibihty that the lensing mass estimates can be subestimated. lOguri et al. 



( I2OO5I ) also investigated the importance of the halo triaxiality and demonstrated that this 
could cause a significant bias in estimating the virial mass and concentration parameter from 
the lensing information. 

In contrast, the concentration parameters obtained either for the single or for the double 
component clusters may be slightly lower ( c. f. Ta bles U an d Bl) , a lthough they agree, within 
the errors, with the expected values . Following IGavazzil (120051 ) . a halo with axis ratio ^ 



1 could expla i n the possi ble small discrep ancies between our model and those inferred by 
Bullock et all J200lh and lEke et all J200lh (c ~ 4 and c ~ 5) for a cluster with 1.1 x 10^^ 
M0 and 1.1 X 10^^ Mq in a f^A-dominated universe. 

A lower concentration changes the mass estimate. At large radii, the effect tends to 
vanish and the weak lensing and true masses are the same (see Figure |8]), but inside the arcs 
(r < 0.1 Mpc) the effect is not negligible producing an overestimate with respect to its actual 
X-ray m ass. But the differe nces in the masses can also be explained as a result of the merging 
process. iRasia et al.l (120061 ) studied a set of five galaxy clusters resolved at high resolution in 



a hydrodynamic simulation, examining the systematics affecting the X-ray mass estimates. 
They showed that for a cluster undergoing a merger, the assumption of hydro-dynamical 
equilibrium led to the underestimation of mass by 20%, and that a (3 model gave even more 
discrepant results with typical deviations of about 40%. In Figure |8] we compare the strong 
lensing mass profile to the X-ray mass profile corrected for the 40% underestimation; we can 
see that with this correction, the X-ray mass estimate agrees with the mass obtained from 
the NFW components at the positions of the arcs, but for smaller radii the X-ray mass is 
significantly lower. It is important to notice that we did not have included the masses of the 
individual cluster components, and thus we are not considering in Figure |8] the total mass 
within the arcs. In addition, the extrapolation of the NFW mass inferred from gravitational 
lensing up to 1 Mpc from the BCG, well beyond the distance of th e arcs, is a factor 2 
lower than the mass shown in Table [1] obtained by iTran et al.l (120051 ). 



NF W ellipsoids have been su ccessfully used for strong lensing modelling in galaxy clus- 



ters (e.g. IComerford et al.ll2006l ). However, the non-linearity of the gravitational lensing 



effect, as well as its high dependence on the core densities, the asymmetries of the mass 
distribution, and the cluster neighborhoods, put hmits to the use of analytical solutions in 



- 15 - 



the ca se of non relaxed systems, such as MS 2053. In this sense, the work of iTorri et al. 



( 120041 ). who have computed the strong lensing effect in a merging cluster with complex and 
irregular mass distribution, have shown that the lensing cross-sections can grow by one order 
of magnitude and that during the merger, the shape of the critical and c austics lines c hange s 
substantially. These topics have been discussed in a recent paper by iFedeli et al.l (120061 ). 
who propose a fast method to calculate cross sections for complex and asymmetric mass 
distributions. 



6. CONCLUSION 

We have modeled the mas s distribution of the cluster MS 2053 using the LENSTOOL 



ray-tracing code developed by iKneibI (119931 ). The fits were performed considering both a 
single and a double NFW profiles, and 23 Galaxy-scale mass components as perturbations 
to the cluster potential. We used ESISs for the individual galaxies and set constraints to 
their parameters using observational data. We measured the arc positions and shapes, and 
the galaxies positions using the HST archive image, and calculated the NFW parameters for 
the single and double galaxy cluster models. Our main results can be summarized as follows: 

1. The XY parameter space of the single ABC model has a bimodal distribution that 
strongly suggests a double cluster undergoing a merger process. The quality of the fits 
also favor a model with a bimodal mass distribution 

2. The models consistently predict a third slightly demagnified counter- image for the AB 
arc system. 

3. We found a candidate for this counter-image near the predicted position. This candi- 
date shows a fuzzy object that can be the image of a distorted galaxy due to the effect 
of gravitational lensing. Besides, it has the same F702W-F814W colors cls circs A and 
B. 

4. Using the strong lensing effect we confirm the asymm etry and the elongation in the 



mass distribution of MS 2053 reported by other authors (IHoekstra et al.ll2002l : iTran et al 



20051 ) . and we estimate the total lensing mass within the arcs to be 4.7 x 10 M^ 



0- 



We find a significant discrepancy between the mass estimates from lensing and X-ray 
measurements, that we attribute to the fact that the cluster is not in hydro-dynamical 
equilibrium. 



5. 



The concentration parameters obtained either for the single or for the double com- 
ponent clusters may be slightly lower than the expected values. These possible small 
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discrepancies can be explained by asymmetries in the mass distribution and the pro- 
jection effects. 

Much more work is necessary in order to understand MS 2053 mass distribution. Future 
spectroscopic follow-up of arcs B and C will provide a test of our model. Also, an extensive 
study of the internal velocity dispersion of the cluster members can be used to further improve 
the accuracy of the mass reconstruction. 
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A. The Einstein angle for a NEW profile 



Consider a spherical NFW density profile acting like a lens. The analytical solutions 
for this lens were given by iBartelmann fll996h and have been studied b y different authors 
dWright fc BrainerdI I2OO0I : IColse fc Kn"eibl I2OO2I : iMeneghetti et al.ll2003[ ). The positions of 
the source and the image are related through the equation: 



= u^ - V(p(u^) = u^ - a{u^) (Al) 

where and are the angular position in the image and in the source planes, respectively. 
a is the defiection angle between the image and the source and (p is the two-dimensional lens 
potential. We introduce the dimensionless radial coordinate x = /u^, where ul = Vg/DoL, 
and Dql the angular diameter distance between the observer and the lens. In the case of an 
axially symmetric lens, the relations become simpler, as the position vector can be replaced 
by its norm. 
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The deflection angle then becomes (iGolse fc Kneibll2002l ): 



a(x) = 4:ks^g{x)ex 



(A2) 



Where g{x) is a function related w ith the surface density inside the dimensionless radius x, 
and is given by (jBartelmann 1996 ): 




(A3) 



Where kg = r.^prritSr^ rrit^ the lensing st rength, is an estimate of the convergence parameter 



(see the discussion in 



Oguri et all 120041 ) . The quantity S 



crit 



{cy47rG){Dos/DoLDLs) is 



the critical surface mass density for lensing, and D^s and Dps are the angular diameter 
distances between the lens and the source, and the observer and the source, respectively. 

If we calculate the deflection angle at radius r^, we obtain x = 1, and equation IA2I can 
be expressed as: 



aix) 



OS 



Hence, the lens equation lAll can be written as: 



(A4) 



u 



u 



1 + 



u 



E.r 



(A5) 



Where u^^ ^ is the Einstein angle that for the case x = 1 is expressed by: 



u 



c'Dos V 2 



(A6) 



For a circularly symmetric lens the Einstein radius is given by = Ana"^ D ls / Dqs ■ There- 
fore, equation IA6I takes the form: 



u 



c'Dos 



(A7) 
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Where = 4 (l + In |) Gr'^pcritSc- Therefore, represents the velocity at radius r^. We 
adopt this characteristic velocity in order to compare the velocity predicted by our models 
with the velocity dispersion measured in dynamical studies. 



REFERENCES 



Bartelmann, M. 1996, A&A, 313, 697 



Bertin, E., & Arnouts, S. 1996, A&A, 117, 393 



Blandford, R.D., & Kochanek, C.S. 1987, ApJ, 321, 658 
Broadhurst, T., Bem'tez, N., Coe, D., et al. 2005, ApJ, 621, 53 
Bullock, J.S., Kolatt, T.S., Sigad, Y. et al. 2001, MNRAS, 321, 559 
Clowe, D., De Lucia, G., & King, L. 2004, MNRAS, 350, 1038 

Comerford, J.M., Meneghetti, M., Bartelmann, M., & Schirmer, M. 2006, A&A, 642, 39 
Covone, G., Kneib, J. -P., Soucail, G. et al. 2005, A&A, submitted astro- ph / 051 1332| 



Delia Ceca, R., Scaramella, R., Gioia, I.M., Rosati, P., Fiore, F., & Squires, G. 2000, A&A, 
353,498 

Dyer, C.C., & Roeder, R. C. 1976, Nature, 260, 764 
Eke, V.R., Cole, S. & Frenk, C.S. 1996, MNRAS, 282, 263 
Eke, V.R., Navarro, J.F., & Frenk, C.S. 1998, ApJ, 503, 569 
Eke, V.R., Navarro, J.F., & Steinmetz, M. 2001, ApJ, 554, 14 
Faber, S.M., & Jackson, R.E. 1976, ApJ, 204, 668 

Fedeli, C, Meneghetti, M., Bartelmann, M., Kolag, K., & Moscardini, L. 2006, A&A, 447, 
419 

Gavazzi, R. 2005, A&A, 443, 793 

Gavazzi, R., Fort, B., Mellier, Y., Pello, R., & Dantel-Fort, M. 2003, A&A, 403, 11 
Ghigna,S., Moore, B., Governato, F., et al. 2000, ApJ, 544, 616 



- 19 - 

Gladders, M.D., Hoekstra, H., Yee, H.K.C., Hall, P.B. & Barrientos, L.F. 2003, ApJ, 593, 
48 

Golse, G., & Kneib, J.-P. 2002, A&A, 390, 821 

Hoekstra, H., Franx, M., Kuijken, K., & van Dokkum, P.G. 2002, MNRAS, 333, 911 

Jing, Y.P., & Suto, Y. 2002, ApJ, 574, 538 

Karoji, H., & Nottale, L. 1976, Nature, 259, 31 

Keeton, C.R., Koclianek, C.S., & Falco, E.E. 1998, ApJ, 509, 561 

Kneib, J.-P. 1993, Ph.D. thesis, Universite Paul Sabatier, Toulouse 

Kneib, J.-P., Ellis, R.S., Small, I., Couch, W.J., & Sharpies, R.M. 1996, ApJ, 471, 643 

Kneib, J.-R, Hudelot, P., Ellis, R.S. et al. 2003, ApJ, 598, 804 

Koopmans, L.V.E., Bolton, A.S., Buries, S., & Moustakas, L.A. 2006, ApJ, 649, 599 
Lacey, C., & Cole, S. 1993, MNRAS, 262, 627 

Le Fevre, O., Hammer, F., Angonin, M., Gioia, I.M., & Luppino, G.A. 1994, ApJ, 422, L5 
Limousin, M., Richard, J., Kneib, J.-P,et al. 2006, A&A, submitted [astro-ph /0612165] 
Luppino, G.A., & Gioia, I.M. 1992, A&A, 265, L9 

Luppino, G.A., Gioia, I.M., Hammer, F., Le Fevre, O., & Annis, J.A. 1999, A&A, 136, 117 
Lynds, R., & Petrosian, V. 1986, BAAS, 18, 1014 

Meneghetti, M., Bartelmann, M., & Moscardini, L. 2003, MNRAS, 340, 105 

Moore, B., Governato, F., Quinn, T. et al. 1998, ApJ, 499, L5 

Narayan, R., Blandford, R., & Nityananda, R. 1984, Nature, 310, 112 

Navarro, J.F., Frenk, C.S., & White, S.D.M. 1997, ApJ, 490, 493 

Oguri, M., Inada, N., Keeton, C.R., et al. 2004, ApJ, 605, 78 

Oguri, M., Takada, M., Umetsu, K., & Broadhurst, T. 2005, ApJ, 632, 841 

Paczyhski, B. 1987, Nature, 325, 572 



-20- 

Rasia, E., Ettori, S., Moscardini, L., Mazzotta, P., et al. 2006, MNRAS, 369, 2013 

Sand, D.J., Treu, T., Ellis, R.S., & Smith, G.R 2005, ApJ, 627, 32 

Sand, D.J., Treu, T., Smith, G.P., & Ellis, R.S. 2004, ApJ, 604, 88 

Smith, G.P., Kncib, J.-R, Ebeling, H.„ Czoskc, O., & Small, I. 2001, ApJ, 552, 493 

Smith, G.R, Kneib, J.-R, Small, I., & Mazzotta, P. 2005, ApJ, 359, 417 

Soucail, G., Kneib, J.-R, & Golse, G. 2004, A&A, 417, L33 

Soucail, G., Fort, B., Mellier, Y., & Pleat, J.R, A&A, 172, L14 

Torri, E., Meneghetti, M., Bartelmann, M., et al. 2004, MNRAS, 349, 476 

Tran, K.V. 2002, Ph.D. thesis, Univ. California, Santa Cruz 

Tran, K.V., Pranx, M., lUingworth, G., Kelson, D.D., & van Dokkum, P. 2003, ApJ, 599, 
865 

Tran, K.V., van Dokkum, P.G., lUingworth, G.D., Kelson, D., Gonzalez, A., & Pranx, M. 
2005, ApJ, 619, 134 

Turner, E.L., Ostriker, J.R, & Gott, J.R., III 1984, ApJ, 284, 1 

Vikhlinin, A., VanSpeybroeck, L., Markevitch, M., Formarn, W.R., & Grego, L. 2002, ApJ, 
578, L107 

Wright, CO., & Brainerd, T.G. 2000, ApJ, 534, 34 

Wuyts, S., van Dokkum, P.G., Kelson, D.D., Pranx, M., & lUingworth G.D. 2004, ApJ, 605, 
677 

Zaritsky, D., & Gonzalez, A.H. 2003, ApJ, 584, 691 



This preprint was prepared with the A AS IM^jX macros v5.2. 



- 21 - 




-30 -20 -10 10 20 



arcsec 



Fig. 1.— HST WFPC2 F814W image of the central region (60" x 60") of the MS2053 cluster, 
North is up, East is left. The arcs A and B are identified, as well as the 23 galaxies considered 
in our models (lower-case letters). The square shows the region enlarged in Fig. [31 
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Fig. 2. — The cluster galaxies and the Faber- Jackson relat ion. The solid line shows the 
correlation for the early- type galaxies of IWuyts et al.l (120041 ) for which the velocity disper- 
sion is available. Dashed line is a similar correlation but for galaxies not included in the 
Wuyts et al. J2004h sample. The tri angle represents galaxy n, whose velocity dispersion has 
been measured by iTran et al.l (120031 ) . 
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Fig. 3. — A detail of the HST observation in the F814W filter of the region shown in Fig. 
[H centered on the counter-image candidate (North is up, East is left). Superimposed on the 
image, the isophot map in arbitrary units shows the shape of a deformed galaxy. 
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Fig. 4. — Two-dimensional distribution of the position of the cluster center for the SINGLE 
ABC model. Two clumps are easily distinguished. We used these clumps to construct the 
DOUBLE ABC model. The points show the results from the different trials of the MCMC 
optimization. 
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Fig. 5. — Projected mass distribution derived from the double component lens model. Note 
the elongation in the cluster mass which traces the distribution of the galaxies. 
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Fig. 6. — HST WFPC2 F814W image of MS2053. The contours correspond to the projected 
surface densities of 2.3, 4.8, 7.3 x 10^° Mq arcsec"^. As well as in Figure [5] we can appreciate 
that the shape of the main component is elongated. 
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Fig. 7. — The projected mass as a function of the aperture radius measured from the BCG for 
the three components of the cluster. The red-doted hne and the green-dashed hne illustrates 
the mass profile for Component 1 and Component 2, respectively. The blue-dot-dashed line 
shows the masses of the individual cluster components. The black-solid line is the total mass 
of the cluster. 



- 28 - 




Fig. 8. — The projected mass for the NFW profile as a function of the aperture radius 
measured from the BCG. The black-sohd hne shows the mass profile derived from the sum 
of Component 1 and Component 2 in the double component model. The vertical blue line 
indicates the radius of the AB arc system. The mass estimate inside this radius is reliable. 
The mass outside the AB arc system has been extrapo lated from the NFW model. The 
dashed line illustrates a (3 model derived from data of IVikhlinin et al.l ( 120021 ). The dot- 
dashed line shows th e same 13 model c orrected for the lack of hydrodynamical equilibrium 
in a merging cluster (iRasia et al.ll2006l ). 
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Fig. 9. — The external critic (continuous) and caustic (dashed) hues for the DOUBLE ABC 
model and the position of the lensed galaxy (shown as a small circle) in the source plane. 
The eUipses are the images predicted by our best model. 
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Table 1. Cluster parameters 



Group 


z 


{Mq) 


(km s^^) 


(km s^^) 




Reference 


MS 2053 
MS 2053 - A 
MS 2053 - B 


0.5866 ±0.0011 
0.5840 ± 0.0005 
0.5982 ± 0.0003 


3.7 X 10^5 
1.1 X 10^5 
1.1 X lO^-* 


1523 ± 95 
865 ± 71 
282 ± 51 


88611^1 


2.9 
3.8 
5.4 


1,2 
1 
1 



''Weak lensing data 
Calculated from lEke et~all 1 I2OOJ) 

References. — (1) iTran et alJ l|2005h : (2) iHoekstra et al.l l|2002h 
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Table 2. Data of the galaxies 



ID Name a 5 F&IAW 0o z ctq Reference 





(J2000) 




(J2000) 




(mag) 






n 






(km s ■"-) 




a 


20 


56 


21 


4 


-04 


37 : 


50 


19.13 ± 0.01 





259 


63.9 


0.583 


292.0 ± 10.0 


1 


b 


20 


56 


21 


5 


-04 


37 : 


54 


20.24 ±0.01 





122 


81.8 





5894 


234.0 ± 23.0 


1 


c 


20 


56 


21 


5 


-04 


38 : 


00 


20.93 ± 0.02 





577 


-41.1 





5724 


151.0 ± 17.0 


1 


d 


20 


56 


20 


1 


-04 


37 : 


35 


20.29 ± 0.01 





368 


- 7.5 





5846 


169.0 ± 14.0 


1 


e 


20 


56 


20 


5 


-04 


38 : 


02 


20.84 ± 0.01 





473 


32.9 





5737 


129.0 ± 15.0 


1 


f 


20 


■ 56 


22 





-04 


38 : 


13 


21.47 ±0.02 





278 


80.2 





5800 


138.0 ± 13.0 


1 


9 


20 


■ 56 


23 





-04 


38 ; 


08 


20.19 ± 0.01 





136 


-74.2 





5854 


124.0 ± 13.0 


1 


h 


20 


56 


22 


1 


-04 


38 ; 


18 


21.12 ± 0.02 





374 


- 9.8 





5835 


181.0 ± 13.0 


1 


i 


20 


: 56 


21 


3 


-04 


37 : 


57 


21.09 ± 0.02 





299 


7.1 





5854 


< 200.0 


2 


k 


20 


■ 56 


19 


6 


-04 


38 : 


05 


19.98 ± 0.01 





867 


-66.5 





5912 


< 200.0 


2 


I 


20 


56 


19 


6 


-04 


37 ; 


54 


21.38 ± 0.02 





270 


-68.3 





5824 


< 200.0 


2 


m 


20 


56 


20 





-04 


37 : 


49 


20.88 ± 0.01 





361 


-10.1 





5813 


< 100.0 


2 


n 


20 


56 


21 


1 


-04 


37 : 


22 


21.34 ±0.02 





240 


76.1 





5764 


86.0 


3 


o 


20 


56 


21 


9 


-04 


37 : 


25 


21.67 ±0.03 





063 


26.3 





5799 


< 200.0 


2 


pa 


20 


56 


23 


2 


-04 


37 : 


50 


20.83 ±0.01 





203 


11.7 





5990 


< 100.0 


2 




20 


■ 56 


20 


9 


-04 


37 : 


47 


21.74 ±0.03 





335 


- 9.0 





5978 


< 100.0 


2 




20 


■ 56 


20 


3 


-04 


37 : 


50 


20.04 ±0.01 





116 


65.5 





5990 


< 100.0 


2 




20 


56 


20 


7 


-04 


37 : 


32 


21.41 ± 0.02 





855 


84.8 





5993 


< 100.0 


2 


t 


20 


56 


20 


9 


-04 


38 : 


06 


21.90 ±0.03 





093 


- 4.4 





5827 






u 


20 


■ 56 


21 


7 


-04 


37 : 


45 


21.04 ±0.02 





546 


-47.6 





5876 






V 


20 


56 


22 


2 


-04 


37 : 


31 


21.99 ±0.04 





217 


22.3 





5892 






w 


20 


56 


20 


8 


-04 


37 : 


28 


22.57 ±0.06 





910 


68.7 





5939 






X 


20 


56 


20 


5 


-04 


37 ; 


26 


21.90 ±0.04 





798 


-48.9 





5782 







^Galaxies that belong to the second cluster component, according to lTran et al.l ||2005| ) 

2 _t2 

°The ellipticities are defined as e = ^^q^p-, where a and b are respectively the semi-major and semi-minor axis of 
the elliptical shape. 

References. — (1) Wuyts et al. 2004; (2) Tran et al. 2005; (3) Tran et al. 2003 . 



Table 3. Best-fitting parameters 



Model 


X 


Y 


e 


e 


c 


ra 




Xpos 






(") 


{") 




n 




(kpc) 


(km s^-*-) 






SINGLE AB 


-4.1 ±3.3 


6.6 ±4.8 


0.24 ±0.15 


104 ± 37 


2.6 ±0.7 


1259 ± 312 


1251 ± 165 


1.3 


4.6 


SINGLE ABC 


-2.5 ±0.9 


9.7 ±3.0 


0.25 ±0.11 


83 ± 4 


2.4 ±0.6 


1265 ± 290 


1136 ± 100 


0.2 


3.4 


DOUBLE ABC 
















1.2 


1.7 


Component 1 


-1.7± 1.4 


7.2 ±2.5 


0.34 ±0.16 


81± 18 


3.1 ±0.9 


634 ± 159 


746 ±114 






Component 2 


-2.3 ± 1.9 


15.1 ± 1.1 


0.27 ±0.19 


89 ±41 


4.9 ±2.1 


202 ± 50 


392 ± 133 







